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HOMENCLATURE 


a. ■ design variables for the ttpper surface, 

^ 1-1,7 

b. - design variables for the lower surface, 

^ 1-1,7 

Cj - section drag coefficient 

- section lift coefficient 

C - section pltchlng~moment coefficient, 
referenced to quarter-chord point 

C - pressure coefficient, (p.~p)/q 
P — 

C *- pressure coefficient corresponding to a 
^ local Mach number of 1 

c - chord 

M - Mach number 

p - free-stream static pressure 
p^ ■ local static pressure 
q ■ free stream dynamic pressure 
Re “ Reynolds number 
s - move direction vector 
X - chordwlse distance 
y - vertical distance 
a - angle of attack 

Subscripts 

Inv " Inviscld 
Is - lower surface 
max - maximum 
min - minimum 
us - upper surface 
vise - viscous 
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THE WORLDWIDE FUEL SHORTAGE has prompted new 
Interest In reducing the fuel consumption of 
existing and future aircraft and a logical 
place to begin Is by reducing drag. At tran- 
sonic speeds drag reduction can be accomplished 
by eliminating or weakening rapid recompresslons 
or shock waves over the aircraft surfaces. In 
particular, the wing profile can often be 
redesigned to achieve drag reduction. 

Several methods are available to aid the 
designer In developing advanced airfoil sec- 
tions (e.g., the hodograph method (1,2), A an 
Inverse method applied In the physical plane 
(3), and a combined Inve.se-dlrect method (4). 
However, the hodograph procedure Is complicated 
and Its application requires extensive experi- 
ence In applied mathematics and theoretical 
fluid mechanics. The Inverse method requires 
an a priori knowledge of the desirable form 
of the pressure or velocity distribution, and 
constraints are not readily Imposed. The 
combination inverse-direct method Is complicated 
and requires a designer In the "loop" to 
monitor and enhance the convergence of the 
optimization process to a realistic airfoil 
shape. 

In this report the numerical optimization 
design technique (5-8) has been extended to 
the design of lifting transonic airfoils. The 
numerical optimization design technique uses 
two existing computer programs; an optimization 
program based on the method of feasible direc- 
tions (9) and an aerodynamics analysis program 
based on an iterative solution of the full 
potential equation for transonic flow (10). 

The optimization procedure can be used to 
design airfoil sections for any speed regime 
from low speed through transonic with geometrlr 
or aerodynamic constraints. Several examples 
of the application of the procedure to the 
design of low-drag transonic airfoils are given. 
It should be noted that the results presented 
here are preliminary and are Intended only to 
illustrate the usefulness and simplicity of 
the technique. 

DESIGN METHOD 

The program organization of the numerical 
optimization design technique used in this 


^Numbers in parentheses designate References 
at end of paper. 


study is shown In Fig. 1. The initial airfoil 
required to start each design problem is obtained 
by fitting a polynomial to that region of the 
airfoil to be modified, or to the entire airfoil 
if a complete new profile is desired. The coef- 
ficients of such a polynomial are the design 
variables perturbed by the optimization program 
to achieve the desired design Improvement. If 
a 6th degree polynomial representation of the 
airfoil geometry is used, the coefficients for 
the starting airfoil are obtained by fitting the 
polynomials to seven points on the upper surface 
and to seven points on th& lower surface of a 
known suitable airfoil section. The coeffi- 
cients, along with the Mach number, angle of 
attack, and constraints are the required inputs 
to start the optimization process. 

If the problem being considered is pitching- 
moment minimization, the optimization program 
perturbs the polynomial coefficients one by one, 
returning to the aerodynamics program for evalua- 
tion of the pitching moment after each perturba- 
tion. After all coefficients have been perturbed 
and the direction of change of the pitching 
moment has been noted for each polynomial coef- 
ficient change, the optimization program then 
calculates the partial derivatives (by one-sided 
finite difference) of moment with respect to each 
polynomial coefficient, thus forming the gradient 
of the pitching moment (VC„). The direction in 
which the polynomial coefficients are changed to 
reduce the nose-down pitching moment is 
(the steepest descent direction). The optiriza- 
tlon program then increments the polynomial 
coefficients one to four times in the direction 
indicated by -7C^. The program searches in this 
direction until the moment starts to increase 
(because of nonlinearity in the design space) or 
a constraint is encountered. If either of these 
possibilities occurs, new gradients are calcu- 
lated and a new direction is found that will 
deciease pitching moment without violating any 
constraints. When a minimum value of pitching 
moment is reached with no violated constraints, 
the coordinates of the final airfoil are printed, 
together with the final pressure distribution 
and aerodynamic coefficients. 

A graphic presentation of a hypothetical 
optimization problem using two design variables 
is shown In Fig. 2. Assume that the optimiza- 
tion problem to be attempted is minimization 
of drag coefficient with constraints on lift 
coefficient and airfoil thickness, and with 
the airfoil geometry described by the poly- 
nomials shown in the figure. In this case. 


the optimization program Is allowed to pertnrb 
only ai and 02 to achieve the desired drag 
reduction. This permits only a limited 
reshaping of the upper surface. Note that the 
thickness constraint Is a straight line (that 
Is, a linear function of the polynomial coef- 
ficients) whereas the lift constraint is non- 
linear. The optimization program can handle 
both linear and nonlinear constraints effici- 
ently. 

The design space Is divided Into two 
regions, an infeasible region where one or both 
constraints are violated, and a feasible region 
where no constraints are violated. Minimum 
drag in the feasible region Is sought. Assumed 
starting values of a^ and a2 are depicted by 
point A. The gradient of drag (VC^) is 
calculated for this point, giving the direction 
of change in aj^ and 82 as B - - 70 ^. In this 
case, both and 02 must be decreased to 
minimize drag. Movement in direction 1 con- 
tinues until the lift constraint is encountered 
at point B. At this point, the gradient of 
lift ( 7 C^) is required along with VC^ to define 
a new direction s. Now a^ must be Increased 
and &2 must be decreased to minimize drag by 
moving along the lift constraint. The optimum 
Is shown as the point where the line of constant 
drag with least value In the feasible region 
Intersects the lift constraint. 

If the starting values of a^ and &2 are 
given by point C (Fig. 2 ), the problem begins 
In the infeasible region. Now a direction B 
Is determined that moves toward the feasible 
region with a minimum Increase In drag, while 
overcoming the violated lift constraint. Such 
a move requires an fjicrease in both a^ and a2» 

When the feasible region Is reached, n direction 
8 is then determined that will move along the 
lift constraint until the optimum is attained. 

A more complete discussion of airfoil design 
by numerical optimization Is given in ref- 
erence 8 . 

DESIGN RESULTS AND DISCUSSION 

All design problems considered during this 
study consisted of modifying the upper surface 
of an existing airfoil section to achieve 
transonic drag reduction. The equation used 
to describe that portion of the upper surface 
to be modified is the following cubic equation: 

y ■ ai(x/c) ^ + k(x/c) + a2(x/c)^ + a3(x/c)^ ( 1 ) 


The coefficient, k, Is used to match ordinates 
at the chordwlse station where the forward 
region of the airfoil to be modified ends and 
the fixed aft region begins. Ulien the entire 
upper surface is to be modified, the coefficient, 

k, is used to fix the tialllng edge bluntness. 

The design variables arc the coefficients, 

®1“®3» exponent, a^. The starting 

value of a/, uhb 0.5 for all initial airfoil 
sections considered. The final value of a^^ 
depends on the amount of leading edge bluntness 
required to achieve drag reduction. Each air- 
foil modification reported here was found to 
need more leading edge bluntncss than exhibited 
by the starting airfoil to achieve the desired 
drag reduction and hence the final value of 
*fas less chan 0.5 In each case. 

While substantial drag reductions were 
achieved by using equation (1) , the develop- 
ment of shock-free sections was difficult 
because of the restricted class of sections 
that can be described by such an equation. 
Additional work is needed to obtain more general 
geometric representations for transonic airfoil 
design. In particular, equations are needed 
that will allow decoupling of various regions 
of the airfoil surface without permitting 
wavlness to develop. The inability to develop 
shock-free sections may not be a serious 
limitation of the method, however, since the 
off-design characteristics of such sections 
are often poor in comparison to weak-shock 
sections. In each case presented here the 
shock strength was sufficiently weakened to 
prevent shock- induced boundary layer separation 
from occurring at the design condition (l.e. , 
the shock-Mach number was reduced to at least 

l. 2 in each case). 

In reference to the figures cited below, 
only upper surface pressure distributions are 
shown when the presentation is clarified by 
omitting lower surface pressures. (There were 
only slight changes in the lower surface pres- 
sures when the upper surface was modified.) 

CASE 1. UPPER SURFACE MODIFICATION OF AN 
APPROXIMATE NACA 6-SERIES PROFIIE - The results 
of modifying the upper surface of an airfoil 
to achieve drag reduction at M ■ 0.7 are shown 
In Fig. 3. The initial airfoil is an approxi- 
mation to a 13-percent thick MACA 6-eerles 
section. The lower surface Is exact and the 
upper surface is represented by equation (1) . 

The drag coefficient given in the figure is 
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obtained by a surface preseure Integral and 
represents the drag due to shock losses in an 
Invlscld flow. Note that the drag coefficient 
is reduced to approximately one-fourth that of 
the initial airfoil, The drag reduction is 
accompanied by a reduction in nose-down pitching 
moment, a small loss in lift, and a redistri- 
bution of the area contained within the profile 
contour. A constraint on the cross-sectional 
area ("volume") of the profile prevented 
further thinning of the section to achieve the 
drag reduction. Note the substantial reduction 
in shock strength and the forward movement of 
the shock position. 

The aerodynamic characteristics of the two 
airfoils of Fig. 3, corrected for viscosity, 
are shown at Reynolds numbers of 1.3 x 10^ 
and 20 x 10^ in Figs. 4(a) and 4(b), respec- 
tively. The boundary layer characteristics 
were calculated by a theory described in ref- 
erence 10. This theory uses an invlscld pres- 
sure distribution as input to the von Kerman 
momentum equation which, when integrated, 
gives the boundary layer displacement thickness. 
This displacement thickness Is smoothed and 
added to the airfoil contour to account for a 
turbulent boundary layer. The displacement 
thickness correction is computed iteratively 
with the flow calculation. As expected, the 
lift and pitching moment coefficients show 
the typical decamberlng effect of the boundary 
layer; namely, the coefficients for viscous 
flow are smaller than the corresponding coef- 
ficients for invlscld flow. Note that the 
drag reduction due to modification of the 
airfoil contour is smaller when the boundary 
layer correction is Included In the calculation; 
nevertheless, a substantial reduction still is 
achieved. The aerodynamic coefficients shown 
in Fig. 4(b) at a Reynolds number of 20 x 10^ 
are somewhat closer to the corresponding 
invlscld values because of the smaller boundary 
layer thickness at the higher Reynolds number. 
Again, a substantial drag reduction is noted. 

A graph of Cj vs Mach number at a ■ 0“ 
for the initial and final profiles of Fig. 3 
is shown in Fig. 5 for a Reycolds number of 
20 X 10®. The drag divergence Mach number is 
Increased by approximately 0.04 by the contour 
reshaping shown In Fig. 3. It Is Interesting 
to note that this type of re-contouring to 
reduce shock drag at transonic speeds does not 
produce a drag penalty for subsonic flow or a 
morn rapid drag rise above the design Mach 
number (M » 0,7). 


Plot* of low speed (M * 0<1) pressure 
distributions for a * 10* are given in Fig. 6 
for both of the airfoils shown in Fig. 3« The 
low speed aerodynamic characteristics were 
evaluated for the airfoil shown In Fig. 6 and 
for the next two profile modifications (Figs. 7 
and 11) to investigate the compatibility of 
section re-contouring for transonic drag reduc- 
tion with low speed, high lift requirements. 

Note that the pressure peak near the leading 
edge on the upper surface Is reduced by more 
than a factor of 2 by increasing the leading 
edge bluntness of the ’initial’ airfoil shown 
in Fig. 3. Moreover, the adverse pressure 
gradient following the leading edge pressure 
peak is substantially reduced, thereby improving 
maximum lift and low-speed handling. 

CASE II. FORWARD UPPER SURFACE MODIFICA- 
TION OF AN APPROXIMATE NACA 6-SERlES PROFILE - 
The results of modifying the forward 35 percent 
of the upper surface of a 12-percent thick NACA 
b-serles section to reduce shock drag at M" 0.72 
are shown in Fig. 7. Again, a substantial 
reduction in drag is achieved along with a 
smaller nose-dovm pitching moment and little 
change in lift coefficient. The slight irregu- 
larity in the pressure distribution near the 
38-percent chord station is due to a slight 
mismatch in the first derivative of the contour 
at the point where the forward section described 
by equation (1) joined the fixed aft section 
of the airfoil. Such irregularities can arise 
because the only boundary conditions imposed 
on the geometry of the upper surface were 
matched ordinates at the end of the forward 
section. 

The aerodynamic characteristics of the air- 
foil sections of Fig. 7, corrected for viscosity, 
are shown in Figs. 8(a) and 8(b) for Reynolds 
numbers of 1.3 x JO^ and 20 x 10^, respectively. 
As noted for the preceding case (Figs, 3 and 4), 
the lift and pitching moment coefficients are 
reduced due to the decamberlng effect of the 
boundary layer. The reduction in drag is 
smaller when viscous effects are included but 
still large enough to warrant consideration of 
such contour modification for retrofit of 
existing aircraft. 

The effect of Mach number on the drag 
characteristics of the airfoils of Fig. 7 is 
shown in Fig. 9 for a Reynolds number of 
20 X 10^. As for the preceding case (Figs. 

3-6), no drag penalty is incurred at subcrltical 


Mach numbers by the type of contour modifications 
considered here. The drag divergence Mach number 
Is Increased by appr /Ximately 0.02 with similar 
drag rise character ..sties above the design point 
(M 0.72) for both airfoil sections. 

The low speed pressure distributions at 
high angles of attack for the two sections of 
Fig. 7 are shown in Fig. 10. As noted in the 
previous case (Figs. 3-6), the Increased 
leading edge bluntness of the modified section 
produced a reduced pressure peak, a smaller 
adverse pressure gradien: . and a rearward shift 
of the pressure peak. All these factors should 
Increase the maximum lift coefficient and 
Improve the low speed handling of aircraft 
incorporating this airfoil modification. 

CASE III. FORWARD UPPER SURFACE MODIFICA- 
TION OF THE NACA 23015 PROFILE - The results 
of a forward upper surface contour modification 
designed to reduce the shock drag of the NACA 
23015 airfoil section at M ■■ 0.7 and a ■ 0* are 
shown in Fig. 11. In this case the Initial 
airfoil shown In the figure and the correspond- 
ing upper surface pressure distribution represent 
the exact NACA 23015 profile rather than an 
approximation obtained by using equation (1), 
as was done for the initial airfoil In the two 
preceding NACA 6-series airfoil modifications. 

In this case the modified section Is almost 
shock- free (C^j ■ 0.0005) and exhibits nearly 
the same pitching moment and lift coefficients 
as the original NACA 23015 section. The desired 
drag reduction was achieved with a smaller 
contour change than that made In the two pre- 
ceding cases because the Initial shock was 
weaker (compare Figs. 3 and 7 with Fig. 11). 

The aerodynamic characteristics of the 
original and of the modified 23015 sections, 
corrected for viscosity, are shown in Fig. 12 
for a Reynolds number of 10 x 10^. Again, the 
same trend as noted for the two previous NACA 
6-series modifications is noted: namely, the 

lift and pitching moment coefficients are 
reduced, relative to the Invlscid values, due 
to viscous effects, and the percent reduction 
in drag is smaller. 

The effect of Mach niimber on the drag 
characteristics of the modified and original 
NACA 23015 sections is shown in Fig. 13 for a 
Reynolds number of 20 x 10^. Tne improved 
drag rise Mach number and drag rise character- 
istics exhibited by the modified 23015 section 
are similar to those noted earlier for the 


modified 6~eerles sections (Figs. 5 and 9). 
However, the subsonic drag of the modified 
23015 section Is predicted to be somewhat less 
than that of the original NACA 23015 profile. 

The low speed, high lift characteristics 
of the modified 23015 profile require a more 
careful analysis than the modified &~serles 
sections shown previously because the original 
NACA 23015 airfoil section exhibits one of 
the highest maximum lift coefficients of all 
NACA profiles that have similar design lift 
coefficients (11). The type of profile modifi- 
cation used In this study can generally be 
relied upon to improve the high lift ciiaracter- 
Istics of 6-serles airfoils at low speed because 
most 6-serlos profiles have relatively small 
leading edge radii. However, when such modifi- 
cations are applied to NACA 4- and 5- digit 
sections to reduce drag at transonic speeds, 
the low speed, maximum lift characteristics 
of the modified profiles may incur a slight 
penalty. Such a penalty can often be minimized 
or avoided, If the high lift characteristics 
are considered during the optimization process, 
by careful application of constraints on the 
design. 

The invlsuiu, low speed pressure distri- 
butions for the NACA 23015 section and modified 
section are shown In Figs. 14 anc 15, respec- 
tively, at four angles of attack from 0" to 15". 
Analysis of these pressure distributions shows 
two interesting facts: first, the modified 

profile exhibits greater suction pressures 
than the original NACA 23015 section at all 
angles of attack shown, except at 0"; and 
second, the rate of growth of C with angle 

"min 

of attack Is smaller for the modified section 
than for the original section. This latter 
fact Is illustrated more clearly in Fig. 16, 
which shows a graph of the ratio of minimum 
pressure coefficient for the modified section 
to minimum pressure coefficient for the original 
section vs angle of attack. Note, that the 
ratio decreases above a ■ 7®. Since the 
stalling angle for the NACA 23015 section is 
between 15" and 18®, depending on Reynolds 
number (11) , it is possible that the pressure 
coefficient ratio may reach a value of 1.0 
before stall, and hence the maximum lift 
coefficient for the two profiles may be similar. 
A wind tunnel test would be the only means of 
providing positive evaluation of the high 
lift characteristics of the two profiles. 


CASE IV. FORMASD UPPER SURFACE MODIFICA- 
TION OF THE NACA 64A109 PROFILE - Th« final 
dras minlnlzatlon problan conaldarad In this 
study modification of tha forward upper 
surface of the NACA 64A109 section for M ■ 0.82 
(Fig. 17). The results of this modification are 
similar to those shovm for the 6-serles sections 
shovm In Figs. 3 and 7; that Is, tha drag Is 
substantially reduced along with a reduction In 
nose-down pitching moment and little change In 
lift coefficient. The off-deslgn character- 
istics and low speed, high lift characteristics 
were not evaluated because of time limitations 
and the fact that the trends would be expected 
to be similar to those shown for the previous 
6-serles modifications. The Initial profile and 
corresponding pressure distribution shown in 
the figure are those of the exact NACA 64A109 
profile. 

CONCLUDING REMARKS 

A technique for achieving drag reduction 
of airfoil sections at transonic speed was 
demonstrated. When the method Is applied to 
drag reduction of NACA 6-serles sections, or 
to other profiles with small leading edge 
radii, the reduced drag Is generally accompanied 
by a reduction In nose-down pitching moment, 
little change In lift coefficient for the 
design condition, and Improved low speed, high 
lift characteristics. 

When NACA 5-dlglt sections, or other 
profiles with larger leading edge radii, are 
modified to achieve drag reduction at transonic 
speeds, It was again found that the nose-down 
pitching moment was reduced; the lift coefficient 
remained about the same at the design condition. 
However, a more careful analysis of the low 
speed, high lift characteristics Is required 
and some additional re-contourlng of the profile 
may be required to insure adequate low speed 
aerodynamics. 

Further work Is needed to achieve more 
flexibility in the geometric representation 
used with the numerical optimization technique 
so that a wider range of modifications can be 
achieved. 

Experimental verification of the type of 
modifications developed during this study Is 
required before such designs can be considered 
for retrofit of existing aircraft or used In 
new aircraft. 


Th« off-dcBlgn characterlstlcB of the four 
elrfoll modifications developed during this 
study were good. Furthermore, It Is expected 
that such off-deslgn behavior is typical of the 
type of low drag, transonic airfoil section 
generated by the design technique presented in 
this report. 

In all cases the modifications were devel- 
oped by use of an Invlscld aerodynamic theory 
to reduce computational cost. The aerodynamic 
characteristics of the initial and final pro- 
files were verified by use of a theory that 
Includes viscosity corrections. 

A typical airfoil modification, Including 
re-contourlng by numerical optimization, 
verification of design by use of a theory 
corrected for viscosity, and evaluation of 
low speed, high lift characteristics required 
approximately 20 min of CPU time on a CDC 7600 
computer. 
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Fig. 7 - Inviscid drag minimization; M * 0.72 
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9 Hach number effect on the drag 
characteristics of the airfoils shown 
figure 7; Re “ 20*10®, a • 0* 



INITIAL AIRFOIL 1.A5 0.0031 -0.0732 



Fig. 10 - Low speed pressure distributions 
for the airfoils shown In figure 7; M * 0.1 
Re = 10x10®, a » 10* 



Fig. 11 “ Invlscid drag alnlmization; M ■ 0.7 
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12 ^ Aerodynamic characteristics of the 
^l^rfoils shown in figure 11 with boundary* 
layer correction; M = 0,7, Re 









Fig. 16 - Ratio of Cp . for the NACA 23015 
airfoil to for the modified 23015 

airfoil: M » 0. 1 
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Fie. 17 - Inviscid drag minimization; M * 0.82 
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